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摘要 : 中 国 新 一 代 频 电 频 谱 日 像 仪 - 明 安 图 射电 频谱 日 像 仪 (Mingantu Spectral 
Radioheliograph, MUSER) 以 高 时 间 、 高 空间 、 高 频率 分 辩 率 工作 在 0.4GHz~15GHz， 为 
太阳 爆发 活动 初始 能 量 释放 区 的 物理 过 程 、 太 阳 电 子 加 速 等 研究 开辟 了 新 的 窗口 。 高 性 能 
高 质量 太阳 成 像 算法 是 MUSER 数据 处 理 流水 线 至 关 重 要 的 研究 内 容 。 参 考 法 国 墨 东 天 文 
台 太阳 干涉 阵 的 数据 处 理 方 法 ， 系 统 讨论 分 析 了 多 尺度 (Multi-Scale) CLEAN 算法 ， 给 
了 适用 于 MUSER 的 多 尺度 CLEAN 算法 参数 ， 并 重点 讨论 了 算法 的 GPU 并 行 。 实 验 结 
表明 ， 改 进 的 多 尺度 CLEAN 在 算法 效率 上 比 基 于 GPU 实现 的 Högbom CLEAN 提高 了 近 3 
倍 ， 有 效 提高 了 整个 MUSER 数据 处 理 流 水 线 的 性 能 。 
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射电 观测 是 研究 太阳 剧烈 活动 的 重要 探测 手段 。 建 于 中 国明 安 图 的 螺旋 天 线 阵 列 - 明 安 
图 射电 频谱 日 像 仪 (Mingantu Spectral Radioheliograph, MUSER) 在 400MHz~15GHz 的 
584 个 频 点 上 对 太阳 进行 高 分 辩 率 成 像 ， 其 观测 将 填补 日 本 野 边 山 日 像 仪 (2 个 频 点 ， 
17GHz, 34GHz) ~ H ra pi HAX (5 个 频 点 ，150~450MHz) 、 俄 罗斯 伊 尔 库 茨 克 射 
日 像 仪 (5.7GHz) 等 在 厘米 分 米 波 段 的 空白 区 域 ， 对 太阳 活动 及 对 人 类 影响 的 研究 起 到 习 
要 作用 口 。 
明 安 图 射电 频谱 日 像 仪 分 为 低频 阵 〈MUSER-L 400MHz~2GHz) ) 和 高 频 阵 (MUSER- 
Tl, 2GHz~15GHz) ， 其 高 时 间 、 高 空间 和 高 频率 分 辨 率 给 数据 处 理 系 统 提出 了 较 大 的 挑战 。 
高 频 阵 和 低频 阵 的 数据 接收 机 每 3ms 分 别 接收 一 次 包含 16 通道 的 观测 数据 。 在 后 续 的 数据 
处 理 过 程 中 ， 如 何 高 效 地 获取 高 质量 的 太阳 图 像 成 为 数据 处 理 系统 的 关键 问题 。 
射电 干涉 阵 成 图 的 主要 过 程 包 含 加 权 (Weighting) 、 网 格 化 〈Gridding) 、 傅 里 叶 变 
换 ， 退 卷 积 等 主要 步骤 ， 其 中 退 卷 积 是 对 脏 图 进行 洁 化 的 过 程 ， 是 成 像 过 程 中 最 耗 时 的 部 
分 。 对 日 像 仪 成 像 来 说 ， 退 卷 积 过 程 也 是 整个 管线 设计 中 最 令 人 关注 的 问题 。 本 文 在 细致 
调研 相关 工作 的 基础 上 ， 进 一 步 对 射电 干涉 阵 成 像 中 的 多 尺度 退 卷 积 方法 进行 了 讨论 ， 给 
出 了 相应 的 实现 方法 。 


1， 相 关 研 究 工作 


射电 干涉 成 像 的 基本 原理 是 天 空 实际 亮度 是 可 见 度 函 数 的 傅 里 叶 变 换 。 由 于 在 实际 观 
测 过 程 中 UV 覆盖 的 不 完全 性 ， 采 样 函 数 是 离散 的 。 自 1974 年 H5gbom 提出 CLEAN 算法 
的 思想 后 ，CLEAN 算法 得 到 了 极 大 的 发 展 并 在 射电 望远镜 的 图 像 重 建 算法 中 得 到 了 广泛 的 
应 用 。 经 典 H6gbom CLEAN 算法 的 关键 是 为 天 空 实际 亮度 提供 一 个 模型 () ， 和 迭代 寻找 亮 
度 的 峰值 位 置 ， 并 如 〈1) 式 进 行 于 加， 其 中 参数 为 和 迭代 过 程 中 的 峰值 亮度 ， 为 峰值 亮度 的 
位 置 ， 由 残留 图 像 的 峰值 亮度 确定 3。 (2) 式 中 的 为 脏 图 亮度 ， 为 脏 束 ， 为 公式 1) 计 
算 的 出 的 到 加 后 的 亮度 值 ， 为 残 图 的 亮度 ， 标 记 迭 代 次 数 。 
， (1) 
. (2) 
Hogbom CLEAN 算法 的 一 些 变 种 算法 ， 如 Clark, Cotton-Schwab i. He KIA IE 
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MEM (the Maximum pee Method, MEM) 随后 被 提出 。Steer Clean 算法 被 用 于 日 本 野 
边 山 日 像 仪 的 成 像 ， 其 洁 化 效率 是 Hgbom 算法 的 近 10 

外。 然而 ， 基 本 的 Högbom CLEAN 算法 由 于 在 处 理 点 源 和 展 源 上 表现 的 优势 ， 一 直 被 
广泛 使 用 。 

但 随 着 新 一 代 射 电 干涉 仪 的 不 断 发 展 ， 在 确保 基本 CLEAN 算法 的 简单 性 和 信 噪 比 的 
前 提 下 ， 多 尺度 成 为 提高 效率 的 方法 之 一 。 Multi-Resolution CLEAN, Multi-Scale Maximum 
Entropy, Wavelets CLEAN 等 多 尺度 算法 有 效 提 高 了 洁 化 效率 ， 但 存在 尺度 大 小 固定 和 计算 
代价 高 等 问题 。 

Cornwell 于 2008 年 提出 的 Multi-Scale CLEANB] 算 法 被 认为 算法 直接 且 具 有 稳定 性 和 收 
MIEGS HFA. Multi-Scale CLEAN 为 减少 算法 的 迭代 次 数 ， 提 高 算法 的 稳定 性 ， 采 用 多 个 
不 同 尺度 大 小 的 卷 积 核 与 脏 图 进行 卷 积 ， 得 到 多 个 卷 积 图 像 ， 在 多 个 图 像 中 寻找 全 局 亮度 
峰值 。Multi-Scale CLEAN 算法 的 基本 流程 如 下 : 

(1) ， 选取 不 同 的 尺度 大 小 ， 与 脏 图 进行 卷 积 ， 得 到 卷 积 后 的 系列 图 像 ; 

(2) ， 寻找 全 局 亮度 最 大 点 ， 记 录 亮 度 最 大 点 的 位 置 和 所 在 的 尺度 大 小 ; 

G) “计算 步骤 (2) 中 的 尺度 与 脏 束 的 卷 积 ， 再 乘 以 增益 值 ， 存 储 计 算 结 果 及 尺度 

大 小 ; 

(4) MM (1) 中 所 有 的 卷 积 图 像 中 减 去 (3) 的 计算 结果 ; 

(5) ”重复 以 上 4 个 步骤 直至 到 达 循 环 冰 值 或 达到 预 设 循环 次 数 ; 

(6) 4% (3) 中 存储 的 结果 与 洁 束 卷 积 并 加 上 参与 图 得 到 最 终 的 洁 图 。 

为 了 控制 算法 的 迭代 次 数 ， 提 高 图 片 质 量 ，Multi-Scale CLEAN 算法 需要 指定 不 同 尺 度 
的 卷 积 图 像 不 同 的 权 值 ， 每 次 迭代 循环 结束 后 ， 残 图 需要 乘 以 所 在 尺度 的 权 值 来 增加 或 减 
少 亮度 幅 值 。 经 过 系列 实验 ，“〈3) 式 中 权 值 被 认为 是 较 合适 的 参数 值 局 。 

， (3) 
其 中 ， 为 权 值 ，a 为 当前 尺度 大 小 ; 为 最 大 尺度 大 小 ; 为 权重 参数 。 

为 了 减 小 迭代 算法 对 增益 的 依赖 ，Cornwell 选择 抛物 线 函 数 作 为 Multi-Scale CLEAN 
的 卷 积 核 函数 如 下 : 


(4) 
其 中 ， 为 球面 波 函 数 ; (x, y 为 像素 点 的 位 置 坐标 ) ; AREA). SHR, mo 由 对 应 
Hégbom CLEAN 中 的 Dirac o 


由 于 高 斯 卷 积 核 与 (4) 式 给 出 的 卷 积 核 函数 仅 在 较 高 动态 范围 下 存在 较 小 的 差异 ， 一 
般 来 说 ，Cornwell 采用 不 同 尺 度 的 高 斯 卷 积 核 。 以 洁 束 大 小 为 基准 ， 卷 积 核 大 小 选取 常用 
的 等 比 数列 ， 即 依次 取 洁 束 大 小 的 1、2、4、8 倍 。 

为 进一步 提高 处 理 速度 ， 近 几 年 利用 图 形 处 理 器 技术 加 快 退 卷 积 处 理 成 为 主要 的 方法 ， 
如 基于 图 形 处 理 器 的 网 格 化 算法 使 得 w-projection 的 效率 较 中 央 处 理 器 提高 近 百 倍 母 :项目 
组 在 前 期 基于 图 形 处 理 器 的 -HE | CLEAN 3S 法 较 中 央 处 理 器 干 的 HH | 效率 提高 了 十 
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时 能 够 独立 地 进行 多 个 洁 化 迭代 。 
法 国 默 东 天 文 台 赤 阳 干涉 阵 与 MUSER 一样， 也 是 一 个 去 阳 专 用 的 射电 干涉 阵 。 其 数 
据 处 理 采 用 子 基 手 频率 带 通 滤波 的 多 尺度 洁 化 方法 54。 该 施法 基 手 两 个 箱 化 ， C4》 没有 过 
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SPEARS, ec sc 目 像 仪 的 天 线 分 布 情况 各 最 长 
基线 情况 一 最 终 在 本 交 实 现 的 算法 基本 原理 如 十 : 
< 在 频 域 tr 平面 十 NR SAN arte Beni St sail 
的 半径 上 2 =r Hy py k phe 7 H 
次 将 计算 最 长 基线 反 过 来 计算 相应 的 涨 波 器 带 通 频率 一 
《27 通过 滤波 器 二 将 原来 的 稀 路 UY 分 布 Bw 进一步 拆 分 为 若干 个 不 同 尺 度 的 独 耻 分 布 ， 
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